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We show that the dynamics of disordered charge density waves (CDWs) and spin density waves 
(SDWs) is a collective phenomenon. The very low temperature specific heat relaxation experiments 
are characterized by: (i) "interrupted" ageing (meaning that there is a maximal relaxation time); 
and (ii) a broad power-law spectrum of relaxation times which is the signature of a collective 
phenomenon. We propose a random energy model that can reproduce these two observations and 
from which it is possible to obtain an estimate of the glass cross-over temperature (typically T g ~ 
100 — 200 mK). The broad relaxation time spectrum can also be obtained from the solutions of 
two microscopic models involving randomly distributed solitons. The collective behavior is similar 
to domain growth dynamics in the presence of disorder and can be described by the dynamical 
renormalization group that was proposed recently for the one dimensional random field Ising model 
[D.S. Fisher, P. Le Doussal and C. Monthus, Phys. Rev. Lett. 80, 3539 (1998)]. The typical 
relaxation time scales like r typ ~ to exp (T g /T). The glass cross-over temperature T g related to 
correlations among solitons is equal to the average energy barrier and scales like T g ~ 2x£oA. x is 
the concentration of defects, £o the correlation length of the CDW or SDW and A the charge or 
spin gap. 

I. INTRODUCTION 

In spite of several decades of intensive experimental and theoretical works related to slow relaxation phenomena, 
important questions regarding the nature of the low temperature phase of spin glasses have remained unsolved. Replica 
symmetry breaking JjJ and droplet theory constitute two pictures that were already available in the 80 's. It is a 
debated question to determine which of these two visions of the problem does apply to laboratory experiments (see 
for instance Q). Recent "memory and chaos" experiments Q suggest that a new type of droplet theory is needed but 
there exists unsolved questions that have been the subject of recent works (see for instance HQ])- Other approaches 
have focussed on the description of out-of-equilibrium ageing dynamics in terms of generalized fluctuation-dissipation 
relations § . These ideas were originally developed in the context of spin glass models but have been used recently to 
discuss different systems such as domain growth dynamics or chaotic flows (9) . It has been shown recently in Ref. (!(]] 
that the long time dynamics of domain walls in the one dimensional (ID) random field Ising model (RFIM) could be 
described by a dynamical real space renormalization group (RG). The dynamical RG is similar to the Dasgupta-Ma 
RG (TTJ that was applied to the ID random Ising model in a transverse magnetic field (13] and to the ID random 
antifcrromagnetic Heisenberg model |l3f| . The Dasgupta-Ma RG was also applied recently to many other random spin 
models |l2|-^9|. The general purpose of our article is to use the dynamical RG to describe the long time collective 
dynamics of disordered charge density waves (CDWs) and spin density waves (SDWs). 

Disordered quasi-one dimensional (quasi ID) systems have known a renewed interest recently since the discovery 
of several inorganic low dimensional oxides such as CuGeC-3 pCj|-p3t (a spin-Peierls compound), PbNi 2 V2 0s |24|] 
(a Haldane gap compound) and Y2BaNi05 p5|-[3C| (a Haldane gap compound). In these compounds the magnetic 
sites can be substituted with Zn (a non magnetic ion). These substitutions generate either 3D antiferromagnetism 
at low temperature (in the case of CuGeOa and PbN^V^Os) or strong 3D antiferromagnetic correlations (in the 
case of Y2BaNiOs). In these two types of low dimensional spin models, the non magnetic defects generate spin-1/2 
moments out of the non magnetic ground state, either in the form of solitons in spin-Peierls systems [j3ij-|38[, or in 
the form of "edge" spin-1/2 moments in Haldane gap systems p9|-p2|. The available theoretical description of these 
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systems is based on the Dasgupta-Ma real space RG jy],§l| describing non magnetic defects in ID and quasi-lD 
geometries p(i| , [44| , fl5[ , and appears to be in good agreement with experiments. 

It has been shown experimentally by some of the present authors that CDWs and SDWs present slow relaxation 
phenomena at very low temperature [[46|-[48| under the form of what has been called "interrupted ageing" in the spin 
glass literature p{|. In fact there are two glass cross-over temperatures. One is due to the freezing of the CDW 
domains and is observed at relatively high temperature mainly in dielectric susceptibility experiments and is defined 
in the usual way for glassy systems as a separation of the so-called a- and /3-relaxation processes (see Ref. Q). We 
are interested here in slow relaxation phenomena occurring at lower temperature that were observed in temperature 
relaxation experiments Q-^8| . Slow relaxation in disordered CDWs and SDWs and antiferromagnetism in the spin- 
Peierls compound CuGeC>3 have the common point that in both cases the physics is related to solitons having a slow 
dynamics in CDWs and SDWs (the microscopic time associated to the reversal of the soliton can be deduced from 
experiments and is r ~ 1 s - see Ref. [ffoj-[f"8] ) , and a fast dynamics in CuGe03 (no slow relaxation has been reported 
in ac-susceptibility experiments - see for instance [ [32) - and specific heat relaxation experiments |5lj| ). The large 
value of To in CDWs and SDWs is due to energy barriers associated to the dynamics of individual solitons, not present 
in CuGe03. By comparison the microscopic time To — 10 -12 s in spin glasses corresponds to the reversal time of 
an individual spin. We show in this article from the analysis of quasi ID strong pinning models that the dynamics 
of disordered CDWs and SDWs can be interpreted as a collective dynamics of randomly distributed solitons. This 
dynamics is similar to a domain growth dynamics in the presence of disorder in which a correlation length is 
increasing with time. Larger objects have a slower relaxation because most of the time the energy barrier associated 
to a correlated pair of solitons that are close in space is larger than the energy barriers of the individual solitons. 



It has been shown in Refs. 52-561] from the analysis of a strong pinning model that there exist slow relaxation 
phenomena associated to independent strong pinning impurities in CDWs. It was also shown in Ref. |5(J that the 
explanation based on independent impurities requires an artificially large concentration of impurities (one impurity 
per unit cell). It is therefore a relevant question to reexamine the experimental data and investigate new mechanisms 
responsible for slow relaxation in disordered CDWs and SDWs. More precisely we address the following questions: 

1. In section ||, we reexamine the slow relaxation experiments in CDWs and SDWs already presented in Ref. jh] Q 
and find two features: 

(i) An "interrupted ageing" behavior (all relaxation times are smaller than a maximum relaxation time r max ). 

(ii) A power-law relaxation spectrum signaling the presence of a broad distribution of relaxation times. 

The strong pinning model with independent impurities proposed in Ref. f5(| has been successful to explain (i) 
but it does not explain (ii). In this article we look for a model that is consistent with (i) and (ii). 



2. We propose in section [II a phcnomenological random energy model (REM) that is compatible with the slow 
relaxation experiments. 

3. We show in section |[y] that the experiments can be qualitatively described as being due to a collective behavior 
in a strong pinning model. A dynamical RG is used to describe the coupling between the solitonic deformations 
of the CDW at different impurities at variance with the previous model proposed in Ref. |5(J where the solitonic 
defects are independent from each other. 

4. We show in section |y| that the experiments can also be qualitatively described by the disordered spin-Peierls 
model proposed in Refs. QHQ. 



II. EXPERIMENTS 

The specific heat and heat relaxation experiments (see Refs. [^6|-^8|) have been performed at the CRTBT-Grenoble 
in a dilution cryostat over the typical temperature range 80 mK - 2 K on sample mass of a few hundreds mg. The 
experiments have been performed with a thermal transient technique, the sample being loosely connected to the 
regulated heat sink via a thermal link. This technique enables us to send energy in the sample for variable durations, 
from a pulse of a fraction of seconds up to a long "waiting time" of 24 h or more. For the thermal transient experiments 
reported here, the procedure is the following. Once the sample is in equilibrium with the heat bath at To, one increases 
slightly the sample temperature to To + ATo, (with ATq/Tq < 10%) during a waiting time t w . The energy source 
is switched off at t w and the thermal transient AT(t,t w ) is recorded until the temperature has relaxed to the initial 
temperature Tq. t is the time elapsed since the waiting time t w . The temperature relaxation AT(t,t w ) depends on 
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FIG. 1. Equilibrium relaxation time spectrum of the incommensurate SDW compound (TMTSF^PFg. Eq. (Q) has been 
used to obtain the spectrum of relaxation times from the temperature relaxation signal. The waiting time is long enough so 
that thermal equilibrium has been reached for all temperatures except T = 123 mK. The long-time tail of the spectrum of 
relaxation times is well fitted by a power- law: P oq (log 10 t) = 6x t~ - 7 for T = 500 mK; P eq (log 10 t) = 3 x t~°' s for T = 400 mK; 



Pcq(log 10 t) = 8 x t~ 



for T = 270 mK. 
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FIG. 2. Out-of-equilibrium relaxation time spectrum of the incommensurate SDW compound (TMTSF)2PF6 at the temper- 
ature T = 200 mK. Eq. (Q) has been used to obtain the spectrum of relaxation times from the temperature relaxation signal. 
Thermal equilibrium has been reached for the longest relaxation time (tw = 11 h). The long-time tail of the spectrum of 
relaxation times is well fitted by a power-law: Pt w (log 10 t) — 2.6 x t~ 0,41 for t w — 1 sec; Pt w (log 10 t) = 5 X t -0,41 for t w — 1 min; 
P tw (log 10 t) = 0.8 x t~ 0A1 for t w = 10 min. 



the value of the waiting time (ageing behavior). If t w is sufficiently large, AT(i, t w ) does not depend on t w anymore 
("interrupted ageing" behavior, see p9[). One can start a new run at To with a different t w . During a series of runs at 
different To below 1 K, the sample is never re-heated above 1 - 2 K. Due to the exceptional stability of the cryostat, 
the reference temperature To can be regulated within ±2 x 1CP 4 over several tens of hours. 

We use a standard procedure to deduce a spectrum of relaxation times from the temperature relaxation AT(t, t w ): 

AT(t, t w ) = J P tw (log 10 r)e-*/ T dlog 10 r. (1) 

An approximate expression for the spectrum of relaxation times can be obtained by replacing exp(— t/r) by the 
(^-function f?(log 10 r — log 10 t), which is justified for a broad relaxation time spectrum: 

dAT(t,t w ) 

Pt w logio * ai T~ ■ 2 

We first consider the incommensurate SDW compound (TMTSF)2PF6. The spectrum of relaxation times 
P tm (log 10 1) obtained from Eq. (|^) is shown on Fig. |l] for equilibrium relaxation and on Fig. || for out-of-equilibrium 
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FIG. 3. Equilibrium relaxation time spectrum of the incommensurate SDW compound AsF6 (a) for T = 150 mK and 
t w = 10 min (□), t w = 1 h (■); (b) for T = 240 mK and t w = 10 min (o), t w = 1 h (•). Eq. (g) has been used to obtain 
the spectrum of relaxation times from the temperature relaxation signal. The long-time tail of the spectrum is well fitted by a 
power-law: P cq (log 10 i) = 1.2 x t~° 6 for T = 150 mK; P cq (log 10 t) = 2.8 x t -1 - 2 for T = 240 mK. 



TaS, - T=l 10 mK - Out-of-equilibrium relaxation 



3 



-0.5 



-1.5 



-2 





1 \„ JA _ 








^°ooco o \ 

^"im__ o - c 

TO, 








m 


° • 
o •• 

• 




■\ 


o 

_ • 


* tw=l sec 
tw=5 min 
tw=20 min 


□ 
■ 

o 




• 


tw=2 h 10 min 
i 


• 

1 


i 


1 2 


3 


4 



Log ]0 [t (sec)] 

FIG. 4. Out-of-equilibrium relaxation time spectrum of the incommensurate CDW compound TaS3 at the temperature 
T = 110 mK. Eq. (H) has been used to obtain the spectrum of relaxation times from the temperature relaxation signal. The 
long-time tail of the spectrum is well fitted by a power-law: Pt w (log 10 t) — 0.7 x t~ ' 32 for t w — 1 sec; Pt w (log 10 t) = 2.5 x t~°' S2 
for t w — 5 min. 



relaxation. In both cases the long time tail of the spectrum of relaxation times is well described by a power-law. A 
similar power-law spectrum is obtained for the incommensurate SDW compound AsF 6 (see Fig |3|) and the incommen- 
surate CDW compound TaS3 (see Figs. | and |). 

It is visible on (TMTSF) 2 PF 6 (see Fig. ||) and TaS 3 (see Fig. ||) that the exponent of the power-law relaxation in 
the out-of-equilibrium dynamics is independent on the waiting time. The effect of the maximal relation time is also 
visible on Fig. where the long time relaxation is faster once equilibrium has been reached. 



III. RANDOM ENERGY-LIKE TRAP MODELS OF DISORDERED CDWS 

We want to describe the slow relaxation experiments in CDWs and SDWs discussed in section || by a REM-like 
trap model similar to Refs. J4j||57j (see Fig. ^). The "trap" energies —E a are independent random variables chosen 
in a distribution p{E a ). The model with p(E a ) — Po{E a ) — 1/T g exp (— E a /T g ) solved by Bouchaud and Dean in 
Ref. jl9| is recalled in Appendix We consider here a trap energy distribution in which there is a maximal energy 
barrier E maK : 

Pi(E a ) = JL eXp( 7^; /rg) with < E a < E max . (3) 
T g 1 - exp {-E max /T g ) 
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FIG. 5. Out-of-equilibrium relaxation time spectrum of the incommensurate CDW compound TaS3 at the temperature 
T — 165 mK. Eq. (Q) has been used to obtain the spectrum of relaxation times from the temperature relaxation sig- 
nal. The long-time tail of the spectrum is well fitted by a power-law: Pt w (log 10 1) = 10 3 x t~°' s for t w — 100 msec; 
P tu] (log 10 t) = 2.3 x 10 3 x t' - 8 for tw = 15 sec; P tw (log 10 i) = 4x 10 3 x t~ ' 8 for t w = 1 min; P tw (log 10 t) = 4 x 10 4 x r 096 
for t w = 20 min; P tm (log 10 t) = 11 X 10 5 x t~ U2S for t w = 90 min and t w = 36 h. 



The model with the trap distribution (||) has two properties related to the experiments discussed in section O: 

(i) "Interrupted ageing" behavior due to the presence of a maximal energy barrier. The "complexity" Q 
(T/T g ) In (r crg /T ) introduced in Ref. Q is equal to f2 = E max /T g . 



(ii) A power-law distribution of relaxation times for t q < r n 

IT 1 
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FIG. 6. Schematic representation of the phase space of a trap model. There are no traps having an energy barrier larger 
than -Emax- The maximum relaxation time is defined as r max = to exp (i5 max /T). 



A. Average and typical relaxation times 

The average relaxation time is defined by T av = J rP(\og 1Q r)d log 10 t. For the model with the trap energy distri- 
bution (H) we find 

r av _ T 1 - exp [- (1 - T a /T) (E max /T g )} 

t° T -T g l-e X p(-(E max /T g )) • {) 
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The typical relaxation time corresponds to the maximum of P(log 10 r) and is defined as r typ = exp [((lnr Q ))], where 
r a is the trapping time in trap a: r a = t exp (E a /T). With the trap energy distribution (^) we find 



rtyp 



exp 



T 
T 



exp (E mgx /T g ) - 1 



(6) 



which reduces to r typ ~ r° exp [T g /T] if E max is large compared to T g . The activated behavior of the typical relaxation 
time is in agreement with previous experimental observations j4§| . The variation of the average relaxation time versus 
l/T is shown on Fig. [7| If _E max is weak there is an activated behavior for all values of the temperature, even below 
T g . If £" max increases there is a strong slowing down of the dynamics below T g and the average relaxation time is 
infinite when E max is infinite in which case we recover the behavior considered in Ref. p9[ . The average relaxation 
time plays a central role in the weak ergodicity breaking scenario for spin-glass models (see Appendix |A|) but does 
not play a relevant role in slow relaxation in CDWs and SDWs because it cannot be deduced from experiments in 
these systems. 



o 




FIG. 7. Evolution of the logarithm of the average relaxation time log 10 (r av /to) versus l/T with T g — 1 and different values 



of E n 



B. Relaxation below the glass cross-over temperature T g 



We use the approximations presented in Appendix [A 2| to evaluate the dynamical correlation function Tli(t,t w ) 
associated to the model having a maximal energy barrier E n 



Hi (4,0 



1 



J max • 



f \ T / T 9 



t 



if t w < t < 

^max (7) 



T \ 2 ( t x 1 ~ T/T » 



n l( M.)-l-^j I— J if*<i«,<w (8) 

IIiO, t w ) ~ 1 if t < r max < (9) 

We make two comments: 

(i) The system has thermalized before t w if t w 3> r max . In this case, the correlation function decays exponentially 
with time: U(t,t w ) ~ exp (— i/r max ). Within the approximation used here, this corresponds to TL(t,t w ) ~ 1 if 
t < r max and II(t, t„) ~ if <> r max . 

(ii) If t, t w <C T max , the system has no time to "experience" the existence of the maximal energy barrier -E max . The 
out-of-cquilibrium relaxation is identical to the model discussed in Appendix |A]. 
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Compound 


T (mK) 




T/T a 


T a (mK) 






(TMTSF) 2 PF 6 


200 


0.4 


1.4 


142 






(TMTSF) 2 PF 6 


400 


0.5 


1.5 


267 






(TMTSF) 2 PF 6 


500 


0.7 


1.7 


294 






AsF 6 


150 


0.6 


1.6 


94 






AsF 6 


240 


1.2 


2.2 


109 






TaS 3 


110 


0.3 


1.3 


83 






TaS 3 


165 


0.8 


1.8 


92 





TABLE I. Values of a C x P corresponding to the experiments in section |]J We have indicated the estimated values of T/T g 
and the estimated value of the glass cross-over temperature T g 



C. Relaxation above the glass cross-over temperature T g 

Relaxation above the glass cross-over temperature is directly relevant to experiments. Above the glass cross-over 
temperature the correlation function deduced from the approximations in Appendix ^ are found to be 



)r 



T \T/T g -1 



n 1 (t,t w )~(^jj ifi<^<r max (11) 

IIi(M™)^(yJ ift<r max <^. (12) 

The spectrum of relaxation times is a power-law if i <§; <§; T max : dl\.\/d\at ~ t 1 ~ T t Tg . The exponent 1 — T/T g 
does not depend on t w . This coincides with the experimental behavior discussed in section [n] (see Figs. || and ||). 
The model can be used to deduce an estimate of the glass cross-over temperature T g from the exponents obtained in 
experiments (see section |n|). T g is related to the exponent a Gxp appearing in the experimental power-law relaxation 
spectrum P tm {\og w t) ~ t~ a " xp through the relation T/T g = 1 + a cxp . The values of a exp and the estimations of T, 
deduced from the correspondence between the experiments and the REM-like trap model have been given on Table 
For the three compounds (TMTSF^PFg, AsFg and TaS3, the exponent a exp is increasing with temperature which 
constitutes a qualitative agreement between experiments and the REM-like trap model. The values of T g are of the 
order of T g = 100 -r 200 mK. We should note that the model cannot be used to describe heat pulse relaxation in 
the regime t ^> t w . Namely Fig. || and Fig. || suggest that the exponent a cxp is the same if t t w or t w <C t while 
Eqs. ([To]) - ( |l2"| ) would predict a different power-law in the regimes t *C tu, and t w <C t. This is an indication that the 
exponential trap distribution with a cut-off is not well suited to describe the short time dynamics. 



IV. COLLECTIVE DYNAMICS IN A STRONG PINNING MODEL 



A. The model 

Let us now consider the microscopic model discussed by several authors in Refs. []52]-Q that is used to describe 
the pinning of a disordered CDW: 

^ ' i 

The first term in Eq. ( [l3|) is the elastic energy, the second term is the interchain interaction, and the third term is the 
pining potential. The charge density wave vector is Q and the impurities are at random positions xi along the chain. 
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B. One strong pinning impurity 



Let us recall the solution of the one-impurity model discussed in Refs. j[5 56|. With Vi — 0, ip(x) is the solution of 
the sine Gordon-like equation 

vp d 2 \p(x\ 

w sm ip(x) — V5(x — xi)smip(x) = 0. (14) 

2lT ox 



The low-energy solutions are dipoles made of a superposition of two solitons 
The width of the soliton is 



tanf^Utan^expP^). (15) 



Zo = \^ L - (16) 
V 2nw 

For commensurate impurities having Qx\ = 0, ip is the solution of cos (-0/2) = Vc /V, where the one-impurity 
pinning threshold is given by V c {1) = ^/2wv f /tt. With V > V c (1) there is one unstable solution (ip = 0) and there are 
two stable solutions ip = ±2 arccos (Vc^ /V). The single impurity solutions are degenerate and their total energy is 
Efol = —2(V — Vc 1 ^) 2 /V. With incommensurate impurities (Qx\ ^ 0), ip is the solution of 

^f sia (t) =Vsm(Qx 1 +<P). (17) 



The dipolar low-energy excitations were discussed in Rcf. |56 for independent impurities close to commensurability 
(Qxi -C 1)- Eq. Jl7t > can be solved numerically if Qx\ is not a small parameter. 

Another way of solving the one-impurity model is to consider the soliton profile given by ( |l5| ) as a variational 
solution parametrized by ip. The energy landscape associated to a single impurity is the following: 

eHI(X ) = - + X {A( Xl ) +B( Xl )X a - (16< - C( Xl ))X$ - B( Xl )X§ - (16< - A( Xl )) A 4 } , (18) 

where Xq — tan (■0/4). The coefficients A(xi), B(xi) and C(xi) are given by 

A{xx) = 2Vsm 2 (Qxi) (19) 
B(x x ) =4Vsm(Q Xl ) (20) 
C(xi) =2V[l + 3cos(Qx 1 )}. (21) 



Minimizing Eq. ( |18| ) with respect to Xo leads directly to Eq. (|17]). For commensurate impurities there are two 
degenerate energy minima separated by an energy barrier. For incommensurate impurities the two energy minima 
are not degenerate (there is a metastable state). 



C. Two strong pinning impurities 



The purpose of this section is to replace a cluster made of two impurities by a single effective impurity. We start with 
the case of two impurities that are close to each other (\x2 — X\\ -C £o)- The opposite limit \x 2 — x±\ ^> £o is discussed 
in a straightforward fashion because in this case the two impurities have an independent dynamics. An approximate 
solution for an arbitrary \x 2 — x±\ can be obtained by interpolating between the two limiting cases \ X 2 — x± \ <C £o and 

\X2 ~ Xl\ > £ - 



8 



1. Variational soliton profile 



We look for variational solutions describing two impurities at positions x\ and x 2 > x\ under the form tan (ip(x) /4) = 
F(x), with 

F(x) = tan ( *f) exp ( Jl^A) + tan ( *f) exp (J^A\ , (22 ) 



where i/'i and ^2 are variational parameters that can be determined by minimizing the total energy obtained from 
the soliton profile p5). If the distance between the two impurities is much smaller than Co (\x 2 — x\\ *C Co) the 



two- impurity energy landscape El o l(X ) is a function of the single parameter X given by 



X Q = tan ( % ) + tan [ ^ ] . (23) 



2. Commensurate impurities in the limit \xi — X\\ <C Co 

Let us now consider the situation where the two impurities at positions xi and X2 are commensurate and the 
distance between the two impurities is small compared to the width of the soliton (\x 2 — Xi\ <C Co)- The total energy 
of the two-impurity system takes the form 

4ol(Xa) = (1 ^| )a K<> - ^ + < * 2 } ■ 

f 2") 

The ground state is such that dE^ oi :(Xq)/9X = 0, which leads to 




B^(X*) = --(y-^o) 2 . (25) 

We make two remarks: 

(i) The pinning threshold associated to two commensurate impurities such that \x 2 I Co is ^"i 2 ^ = wCo = , 
equal to two times the pinning threshold associated to a single impurity. The interaction between impurities 

(2) 

increases the pinning threshold and favors the so-called "collective pinning" regime corresponding to V < Ve , 

(2) 

as opposed to the "strong pinning" regime corresponding to V > V c . 

(ii) In the limit of strong impurity pinning (V 3> u?Co); there are two energy minima corresponding to Xq = 

(2) 

tan(?/>/4) + tan(-0'/4) = ±1 and having an energy EI q [(Xq) — —AV. These two minima are separated by the 

(2) 

saddle point X = tan(^>/4) + tan (V> /4) = having an energy El o {.(X' ) — 0. If V > wCo and \x 2 — xi\ <C Co, 
the energy barrier associated to the two impurity system is AE^ = _E to t(^o) — -Etot(^o) = 4V, equal to two 
times the energy barrier associated to a single impurity: AE^ — 2V (see section IV B| ). This shows that the 
interaction between impurities makes the system more glassy. 



3. Incommensurate impurities in the limit \X2 — xi\ >C Co 



The energy landscape associated to two incommensurate impurities at a distance much smaller than the soliton 
width (\x 2 — Xi\ -C Co) can be reduced to an effective single impurity energy landscape: 

(X Q ) = - + * {A + B Xo - (l6<o - <%) X% - B X* - (l6< - A ) * 4 } ■ (26) 

The coefficients Aq, Bq and Cq are obtained as the sum of the coefficients associated to the single impurity energy 
landscapes given by Eqs. © - @: A = A(xi) + A(x 2 ), B Q = B(x x ) + B(x 2 ), and B = C(n) + C(x 2 ). The 
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form ( pq ) of the two-impurity energy landscape can be understood as follows. The elastic and interchain energy of 
the two-impurity model are identical to the one-impurity model if one uses the variables Xq = tan (■0/4) for the 
one-impurity system and Xq = tan(^/4) + tan(^///4) for the two- impurity system. The pinning energy is additive 
since the phase of the CDW at x\ is approximately equal to the phase of the CDW at X2 if \%2 — <C £o- 




FIG. 8. Energy landscape E^l (X) versus X given by Eq. ( |li| ) of two isolated impurities being far apart and such that 
Qx\ — tt/8 (o) and Qx2 = 7r/4 (□). If these two impurities are at a distance much smaller than £o the two-impurity system 
can be represented by an effective impurity having an energy landscape E^l(X) (see Eq. (^6|)). The energy landscape E^ t {X) 
is also shown on the figure (■). 



4-. Limit \xi — K2I 2> £0 



If the two impurities are at a distance much larger than the soliton width (\x2 — x±\ £0) the two impurities 
have an independent dynamics. The two impurities are characterized by the coefficients (A\, B\,Ci) and (A2, B2, C2) 
(see Eq. (|||)). The two decoupled impurities can be replaced by a single impurity characterized by the coefficients 
(Aoo , Boo , Coo ) ■ The coefficients (Aoo, Boo, Coo) correspond either to (A\, B\,C\) if the impurity at position X\ has 
the longest relaxation time or to (A2, B2, C2) if the impurity at position X2 has the longest relaxation time. 

5. Interpolation between the limits |asa — xi\ <C £0 and \xi — xi\ ^> £0 

An arbitrary value of the distance between the impurities can be treated by interpolating between the two limiting 
cases \x2 — x\ | <C £0 & n d \x2 — x\ | 3> £0 discussed above. Namely we suppose that the energy landscape at an arbitrary 



distance \x 2 — X\\ is still given by Eq. (26) 



E^l(X ) = - (1+ ^ 2)2 {A + BX - (l6w£o - C) X 2 - BI 3 - (l6< - i) X 4 } , 



(27) 



and that the coefficients A, B and C interpolate between the solutions already obtained in the limits \x2 — xi \ <C £0 
and \x 2 — x\\ > £ : 

i = ioo + (io-ioo)exp^- |a:2 ~ o Xl1 ^ (28) 
B = Boo + (Bq - Boo) exp ( J X *- X A (29) 



Co 

C = Coo + (Co - Coo) exp (^- |X2 ~ o Xl1 ^ . (30) 
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The correlations mediated by the gaped medium decay exponentially with distance and this is why we use an expo- 
nential interpolation in Eqs. ( p8| ) - (^0|). 

D. Dynamical RG 



(a) 









ED ■ — ED • 


— h H- 





(0 -^--^ — ^ ^3 



FIG. 9. Schematic representation of the formation of correlated objects in the collective dynamics, (a) corresponds to 
decoupled impurities (which is the initial condition of the RG flow at time to)- (b) corresponds to a time ii > to- (c) 
corresponds to a time ti > ti > to- The correlation length is increasing with time: £(£2) > £(ti) > £(£())• 



Distribution of relaxation times (Q=1.41) 




Log 10 ( x/T ) 

FIG. 10. Distribution of the logarithm of the relaxation times with an incommensurate CDW wave vector (Q = y/2~ 1.41). 
The distribution of relaxation times has been obtained with N = 500 impurities and averaged over 100 realizations of disorder. 
The parameters of the strong pinning Hamiltonian given by Eqs. ( |l3| ) are Vf = 1, W = 5 X 10~ 6 , Vi = 56 X 10 -3 for all sites 
i. The correlation length in the simulation is £ = 178ao, with ao the lattice spacing. The one-impurity pinning threshold is 
V c = 5.6 x 10 -3 . The impurity concentration is x — 0.05 <C 1. The temperature in units of vf corresponding to the different 
curves are: T = 2.6 x 10~ 3 (□), T = 1.3 X 10" 3 (■), T = 0.67 x 10~ 3 (o), T = 0.35 x 10" 3 (•). The parameters are such that 
the width of the soliton £ and the ratio T/vf have the correct order of magnitude. For instance £ ~ 4000A in CDW compounds 
and ao ~ 3~7A (a = 3.34A in TaS 3 and a = 7.3A in (TMTSF) 2 PF 6 ). We deduce that £ ~ 500-h 1000a which is compatible 
with the value of ^ used in the simulation. The value of the impurity concentration x is also realistic in the sense that it has 
been shown experimentally that the introduction of 0.5% of extrinsic impurities does not modify the slow relaxation properties 
so that x is presumably larger than 0.5%. 

Now we discuss the collective dynamics of the model defined by Eq. ( |l3[ ) . The method is similar to Ref . jl0| and 
consists in eliminating the fastest degrees of freedom. There is a smallest relaxation time r m i n (being the smallest of 
the relaxation times of individual impurities) that increases in the course of the RG. Assuming a broad distribution of 
relaxation times, the density of relaxation times is given by the number of impurities yo(T m i n )(5r mn i n that are eliminated 
as the smallest relaxation time is increased from T m j n to T m j n + (5T m j n . 

We consider that the initial condition of the dynamics is a quench from high temperature. The initial condition of 
the RG flow corresponds to uncorrelated impurities. The experiments discussed in section [o] correspond to a different 
initial condition but it is expected that the two types of initial conditions can be used to discuss the qualitative 
physics. 
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Distribution of relaxation times (Q=0) 




Log 10 ( T/T ) 

FIG. 11. Distribution of the relaxation times in the presence of commensurate impurities having Q = 0. The Hamilto- 
nian parameters are identical to Fig. |lo[ The temperature in units of vp are: T = 2.6 x 1CF 3 (□), T = 1.9 x 1CF 3 (■), 
T = 1.3 X 1(T 3 (o), T = 0.95 x 10~ 3 (•). 



Independent impurities (Q=1.41) 



T=2.6 10-3 
T=1.3 10-3 
T=0.67 10-3 
T=0.35 10-3 




100 150 200 250 300 

Log ,„( T/T„) 



FIG. 12. The same as Fig. nCn but with non interacting impurities (this is the initial condition of the RG flow). 



The energy landscape of an impurity at site Xi is characterized by the coefficients (A(xi), B(xi), C(xi)) (see Eq. |l8| ) 
and by a relaxation time Tj. We note by r m i n = Minjri} the smallest of these relaxation times corresponding to 
impurity Iq. The impurity io has two neighboring impurities: one at the left (site zl) and one at the right (site in). 
The relaxation times of the impurities at sites %l and in are tl and tr. There are two possibilities to eliminate the 
impurity at site io- 

(i) Transform the two impurities at sites io and ih into an effective impurity at site i' L having a relaxation time t' l . 

(ii) Transform the two impurities at sites io and in into an effective impurity at site i' R having a relaxation time t' r . 

The transformation (i) is implemented if t' l < r' R . The transformation (ii) is implemented if t' r < t' l . 

The distribution of relaxation times is shown on Fig. ^ for an incommensurate CDW wave vector and on Fig. ^ 
for a commensurate CDW wave vector. The two systems are qualitatively similar in the sense that (i) the collective 
dynamics generates a broad spectrum of relaxation times; and (ii) there is a maximum relaxation time r max . At long 
times, the spectrum of relaxation times is approximately a power-law, which is not against the experiments discussed 
in section [n|. For comparison we have shown on Fig. the distribution of relaxation times of non interacting strong 
pinning impurities with the same parameters as Fig. HQ There is already a distribution of relaxation times associated 
to independent impurities which is due to the fact that the energy landscape associated to an impurity at position x\ 
depends on the coordinate x\ through Eqs. ( |l9| ) - (^l|). But the distribution of relaxation times of interacting impurities 
is obviously much broader than the distribution of relaxation times of non interacting impurities (see Figs. O and O). 
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Given that there are evidences that a broad spectrum of relaxation times is present in the slow relaxation experiments 
either for commensurate of incommensurate systems, we conclude that the long time dynamics of disordered CDWs 
and SDWs is a collective phenomenon that can be qualitatively captured by the dynamical RG of the strong pinning 
model. The small time dynamics can be well approximated by independent impurities corresponding to Fig. |l2| on 
the condition that the initial condition is a quench from high temperature. 



V. COLLECTIVE DYNAMICS IN A DISORDERED SPIN-PEIERLS SYSTEM 



A. The model 




-Non magnetic deffects 
Spin- 1/2 moments 



A 




Singlets 



FIG. 13. Schematic representation of non magnetic substitutions in a spin-Peierls system. We assume that the spin-1/2 
solitonic moments have a slow dynamics. 

Now we consider the model introduced in Refs. J30|,fl4|,fl5) to describe non magnetic substitutions in spin-Peierls 
and Haldane gap systems (see Fig. |l3| ). In this model non magnetic impurities in a dimerized system generate 
solitonic spin-1/2 moments out of the non magnetic singlets. Two spin-1/2 moments at distance I are coupled by 
an antiferromagnetic exchange J(l) — Aexp(— £/£o) that decays exponentially with distance. £ is the correlation 
length associated to the gaped background. There is an energy barrier associated to the dynamics of an isolated 
soliton in CDWs and SDWs (sec for instance Fig. ||) which explains that the microscopic time r is of order of 
1 sec (see p6|-^8|). The model proposed for non magnetic substitutions in CuGe03 can thus be "transformed" into 
a model of slow relaxation in CDWs and SDWs just by changing the time scale To associated to the dynamics of 
individual solitonic spin-1/2 moments. From this analogy we deduce an expression of the glass cross-over temperature 
of disordered CDWs and SDWs. 

We represent the solitonic spin-1/2 degrees of freedom in the spin-Peierls system by Ising spins distributed at 
random in ID and use a Glauber dynamics. The Ising spins <ii interact with the Hamiltonian 



where the exchange Jij decays exponentially with distance (see Refs. |p4| , ^5| ): 



Ji j = A exp 



(31) 



(32) 



A is the spin gap, £o is the correlation length associated to the CDW without disorder (see Eq. |l6| for the expression 
of £o m the strong pinning model), and dij is the distance between the Ising spins at sites i and j. Since there is 
no frustration the dynamics of the ferromagnetic and antiferromagnetic models are equivalent and we use here the 
ferromagnetic model. The single spin flip Glauber dynamics of the model defined by Eqs. (31) and (|3^) is given by 
the master equation |5q] 



N 



i{a}P({a},t) 



N 

i=i 



-Oi, a N }P ({oi, 



(33) 
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FIG. 14. Schematic representation of the Ising model. The spins <Ji are distributed at random in ID and interact with the 
Hamiltonian (pit). 



where the transition rates take the form 



Wi{a} 



!1 

2 



1 — cr, tanh I (3Jij <Jj 
jev(i) 



(34) 



where V(i) is the set of neighbors of site i. The form of the transition rates given by Eq. ( |34| ) ensures that the detailed 
balance is verified: PB{a}wi{a} = Ps{(7'}wi{ai, cr}, where Pb{u} is the Boltzmann distribution. 



B. Two-spin model 



Following Refs. |30 44 Ij| we consider the model defined by Eqs. ( |3l| ) and 
distance I. The distribution of distances is P(l) = xexp (— xl) and the exchange (|32|) becomes J{1) 



for two Ising spins o\ and <J2 at 
Aexp(-//e ). 

The relaxation time associated to the two-spin cluster can be obtained from the Glauber matrix defined by (^3|): t(1) — 
To [1 — tanh (J(Z)/T)] . Since we consider the long-time behavior, we use the approximation r(Z) : 
The typical relaxation time follows an Aharenius behavior: 



■exp[2J(Z)/T] 



TO 



= exp ((In 



2x£ A 
1 + 



The activation energy in Eq. ( pq ) is equal to the average energy barrier given by 2((J(Z))) = 2x£o/(l + %£o) 
average relaxation time is given by an Aharenius law in limit of a small dilution of impurities (x^o *C 1): 



III . 

TO To 



4^ 6XP ("T 



(35) 



The 



(36) 



The activation energy in Eq. (|36| ) is equal to the largest energy barrier given by 2Max[J(Z)] = 2A. 

Comparing the typical relaxation time in the REM-like model and the two-spin model (see Eqs. (||) and (36)) we 
obtain an estimate of the glass cross-over temperature in terms of the microscopic parameters: 



T = 



2x£ A 
1 + xJo 



(37) 



which is equal to the average energy barrier. The estimate of T g given by Eq. IfVA) is similar to the estimate of the 



Neel temperature of doped low dimensional oxides discussed in Refs. 
by correlations among the solitons. 



45|. In both cases the physics is controlled 



C. RG of the disordered ID Ising model 



Let us consider two Ising spins 0\ and 02 coupled by a ferromagnetic exchange J and replace these two spins by 
an effective Ising spin a. We note T\ and t-i the relaxation times of the spins o\ and 02 and we note T\ = 1/t± and 
T2 = I/V2 the transition rates. The Glauber matrix of the two-spin system can be diagonalized by forming symmetric 
and antisymmetric combinations of the occupation probabilities (]6l| . We deduce the transition rate f of the effective 
spin a: 

r = i I n + r 2 - yVi - r 2 ) 2 + A ri r 2 tanh 2 ((3 J) 1 . (38) 
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FIG. 15. Representation of the RG transformation. The pair of spins 
having a relaxation rate f — 1/f given by Eq. (Bg). 



i 0+ 2 

(o»a> fio+i) i s replaced by the effective Ising spin <r 



The couplings Jl and Jr between the effective spin a and its neighboring spins (see Fig. |T^) can be obtained by 
equating the partition function of the four-spin system and the partition function of the three-spin system, leading to 



Jr=T arg cosh \[2\J cosh (/? J) cosh (/3Jr) 



JL=Targcosh \/2-\/ cosh (/? J) cosh {(3Jl) 



(39) 
(40) 



The RG transformations can be iterated by eliminating the smallest relaxation time r m j n = Min{r^ ; , t^\^} which 

is either the relaxation time associated to a single spin or the relaxation time t^ +1 associated to the pair of spins 
(i, i + 1). The distribution of logarithm of the relaxation times is shown on Fig. [lfjl Because of the collective dynamics 

Random Ising model 

—r 1 1 1 — I — I — I — r-| 




FIG. 16. Distribution of the logarithm of the relaxation times of the random Ising model defined by ( |3l| ) and (^i|) with 
the parameters A = 1, £ = 10, and x — 0.05, with four different temperatures: T = 1 (□), T — 0.5 (■), T = 0.25 (o) and 
T = 0.125 (•). The "microscopic" time scale is to = 1. 



there is a broad distribution of relaxation times and a maximum relaxation time r max , which is compatible with the 
experiments discussed in section [n|. The disordered spin-Peierls system discussed in this section has thus the same 
behavior as the disordered strong pinning model discussed in section In the case of the Ising model (see Fig. |l^) 
the RG generates only time scales that are larger than tq. In the case of the disordered strong pinning model the RG 
generates time scales that are also smaller than To (being the relaxation time of individual solitons) which explains 
the origin of the maximum in the relaxation time spectrum (see Fig. [To|). 
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VI. CONCLUSION 



To conclude we have discussed the collective dynamics of randomly distributed solitons and shown that this approach 
was relevant to disordered CDWs and SDWs. We have proposed a phenomenological REM-like model that in good 
agreement with experiments. Since the article was already summarized in the Introduction, we end-up with open 
questions: 

(i) We assumed that the solitons are decoupled from each other in the initial condition of the RG. This corresponds 
to a quench from high temperature. The initial condition relevant to the experiments in section || is different 
because the dynamics starts from an equilibrium state and a small temperature variation is applied. It has been 
pointed out that there can be a specific dynamics associated to the temperature dependence of the equilibrium 
correlation length |^32| and it is thus a relevant question to examine similar issues in CDWs and SDWs. 

(ii) We have assumed that only charge or spin degrees of freedom were relevant in CDWs or SDWs. The interplay 
between charge and spin degrees of freedom raises many questions already in pure systems (see for instance p3| ) . 
It is an open question to examine the physics associated to defects in the presence of a charge and a spin sector. 
The interplay between charge and spin degrees of freedom could explain the effect of a magnetic field on the 
slow relaxation properties where it was observed that a weak magnetic field has an effect on the relaxation time 
spectra [6^ ]. 

Finally, it would be interesting to test the existence of slow dynamics effects in other types of experiments, such as 
dielectric response. Indeed, it is suggested in Ref. j5(| that the residual so-called u /3q process" could be the link to 
the low-T, slow heat relaxation phenomena. Dielectric experiments in the temperature range of 1 K and below are 
planned to test this hypothesis. 
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APPENDIX A: WEAK ERGODICITY BREAKING IN A REM-LIKE TRAP MODEL 

In this Appendix we discuss a REM-like trap model having p(E a ) = po(E a ) = Jj- exp (—E a /T g ) that was introduced 



in connection with the so-called "weak ergodicity breaking" property in Ref. j49| . We first recall in section A 1 the 
exact solution used by Bouchaud and Dean in Ref. |^9| in their discussion of the weak ergodicity breaking property. 
We give in section |A 2| a set of approximations that can be used to recover the weak ergodicity breaking property. 



Similar approximations are used in the main body of the article (in section III) for the REM-like trap model of 
disordered CDWs. For this model relevant to disordered CDWs it is not possible to use the same exact solution as 
Bouchaud and Dean and this is why we are lead to use a set of well controlled approximations. 



1. Weak ergodicity breaking via an exact solution 

The authors of Ref. |l9| could demonstrate the so-called "weak ergodicity breaking" property in a model having 
p(E a ) — po(E a ) — ijr exp (—E a /T g ). We note P a (t w ) the probability to find the system in trap a at time t w . The 
evolution of the system is given by the Glauber dynamics Q 

. N N 

-j—Pa(t w ) = - ^ W a ^pP a (t w ) + ^ Wp-taPp^), (Al) 
w p=l,P^a /3=l,/3^a 

where the transition rates are given by w a —>@ — exp (— (3E a ) and wp^, a — roexp(—f3Efj), with tq — 1/ro the 
microscopic time scale and (3 — 1/T the inverse temperature. To solve Eq. flAl| ), it is convenient to make a Laplace 
transform with respect to the waiting time: 



p + OO 

P a (E)= Edt w exp(-Et w )P a (t v 
Jo 



1G 



In Laplace transform, the Glauber dynamics equation (Al) becomes 

EP a {E) + Nr e-^P a (E) = ^ + £ r e-^P p (E). 
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(A2) 



The solution of Eq. (A2) is found to be 



P*(E) 



■fe( r «) 
Efl/s( r /3) : 



(A3) 



with /b(t q ) = Er a /(1 + Er a ), and with r a = exp (/3E a )/(Nro) the average trapping time in trap a. The expression 
of the Laplace transform of the dynamical correlation function is given by 



U (t,E)= PQ (E a )P a (E)exp[-t/T a ]dE a . (A4) 
Jo 

The presence of a dynamical glass transition can be seen from the divergence of the average relaxation time 

— = «-» = / ME a )-dE a = — — (A5) 

TO T Jo T 1 - 1 a 



while the typical relaxation time 



r typ 



TO 



exp 



exp 



T 



follows an Aharenius behavior. Using (A3), we obtain the correlation function 

C°°P0{T)fE{T)e-^dT 



' Po(j)f E (T)d7 



(A6) 



(A7) 



The denominator can be evaluated by replacing /e(t) by /o(t) defined as /o(t) = Et if r < 1/E and /o(t) = 1 if 
t > 1/E. This leads to 



+ OG 



Po(T)fE(T)dr 



-TIT, 



_i_{a (£T0) ^._ £T0 } 



"TO ^---1^9 

Since the waiting time is large compared to to, one has Etq *C 1 and therefore 



Po(t)/e(t)c!t 



T(l - T/T 9 



To evaluate the inverse Laplace transform of Eq. (A7), Bouchaud and Dean used the change of variable u = Je{t), 
from what they deduced the exact result 



IIo(i, t w ) 



1 



T 



{l-uf/^u-^du, 



(A8) 



which is valid if T < T g . In the limiting cases t/(t + t w ) <C 1 and t/(t + t w ) ~ 1, Eq. ( A7) reduces to 



n (t,t w ) ~ i 



T ( t 



n (t,t w )~ 1- — 



■*5 \ ^ H~ ^ 



T 



l-T/Tg 



if t/(t + t w ) < 1 



T/Tg 



if t/(t + t W ) ~ 1. 



(A9) 
(A10) 



The correlation function IIo(£, tu,) tends to zero if f^, is finite and t — * +oo while it tends to unity if t if finite and 
t w — > +oo, which constitutes the so-called "weak ergodicity breaking" property Ea]. 
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2. Weak ergodicity breaking via an approximate solution 



Now let us give a set of approximations that can be used to recover the weak ergodicity breaking property. These 
approximations will be applied to another model in section II] in a situation where we cannot use the change of 
variable leading to (A8). The approximations are the following: (i) We replace /e(t) by fo(r) in the expression ( A7) 
of IIo (t, E); (ii) We replace the exponential exp {—t/r) by the ^-function 9(r — t). We deduce the Laplace transform 
of the dynamical correlation function 



T 



Uo(t,E) ~ 1 - —{Et) l - T / T ° \iEt<\ 
tl (t,E) ~ (1 - L)(Et)- T ^ if Et > 1. 

■L a 



(All) 
(A12) 



The inverse Laplace transform is evaluated within the same approximations. Namely we replace exp (— t/r a ) by 
9(r a — t). This leads to 



U (t,E) 



l/E 



Edt w Ilo[t, t v 



from what we deduce 



LTo(i, t v 







d(l/E) 

The final form of the correlation function is found to be 



E 



n (t,t w ) 

LTo(i, t w ) 



1 



T \ 2 ( t x 1 " T/T f 



1 1, 



E=l/t u 



1 'E 



— J if tw < t 



(A13) 
(A14) 



The approximate correlation functions given by (A13) - (A14) reproduce well the asymptotic behavior of the exact 
solution given by ( A£ ) - ( A10 ) except for the prefactors that are not relevant to our discussion. The qualitative physics 
of the model (namely the weak ergodicity breaking property) can thus be reproduced from these approximations. We 
use the same approximations in the main body of the article for the model having the trap energy distribution (|^) 
for which we cannot use the exact solution anymore. 
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